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ABSTRACT 


Three aspects of barium ion cloud dynamics are discussed. First, 
the effect of the ratio of ion cloud conductivity to background ionospheric 
conductivity on the motion of barium ion clouds is investigated and com- 
pared with observations of barium ion clouds. This study has led to the 
suggestion that the conjugate ionosphere participates in the dynamics of 
barium ion clouds. Second, analytic work on the deformation of ion clouds 
is presented and found to be qualitatively different from that found recently 
by Simon and Sleeper. Third, the linearized stability theory of Linson and 
Workman has been extended to include the effect of the finite extent of an 
ion cloud as well as the effect of the ratio of ion cloud to ionospheric con- 
ductivities. 

The stability properties of a plasma with contra -streaming ion 
beams parallel to a magnetic field are investigated. The results are in- 
terpreted in terms of parameters appropriate for collisionless shock waves. 
It is found that this particular instability can be operative only if the up- 
stream Alfven Mach number exceeds 5. 5. 
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I. INTRODUCTION 


During the past year significant progress has been achieved in both 
areas of study under the present contract. The principal results obtained 
are reviewed in the following two section. The work reported in Section II 
on the dynamics of barium ion clouds was performed by L. M. Linson 
while the work reported in Section III on streaming instabilities was 
performed by K. Golden, L. M. Linson, and A. Sivasubramanian. 

As a result of Project Secede sponsored by the Advanced Research 
Projects Agency during which barium ion clouds formed from different 
size payloads were released at various altitudes, it was possible to test 
theoretical ideas concerning the motion of barium ion clouds. This study 
has led to the suggestion that the conjugate ionosphere participates in the 
dynamics of barium ion clouds as discussed in Section II. A and presented 
at the Fall Meeting of the AGU in San Francisco, California, in December 
1971. Analytic work on the deformation of ion clouds is discussed in 
Section II. B where it is shown that a recent treatment by A. Simon of 
the University of Rochester is incorrect. A note to that effect will be 
submitted shortly to the J. G. R. for publication. Unique observations of 

ion cloud striations made during Project Secede II by looking up the 
magnetic field have led to an extension of the linear theory presented 
earlier by Linson and Workman. This analysis is given in Section III. C 
and was presented at the Annual Meeting of the Division of Plasma Physics 
in Monterey, California, in November 1972. 



In Section III we report the results of a study of the stability properties 
of a plasma with contra -streaming ion beams and interpret these results in 
terms of parameters appropriate for collisionless shock waves. We find 
that this particular instability can be operative only if the upstream Alfven 
Mach number exceeds 5. 5. A paper based on this work is being prepared 
for submission to the Physics of Fluids. 
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II. BARIUM CLOUD STUDIES 


A. INTRODUCTION 

We have previously^ described some of the results that can be 

obtained from a simple two-dimensional model for barium ion clouds 

2 

released in the ionosphere. Recently, Zabusky, Perkins, and Doles 
have given a careful derivation of the coupled nonlinear differential 
equations that have been useful for describing the dynamics of barium ion 
clouds. These equations can also be obtained by starting with simplified 
equations of motion for electrons and ions in the limit of zero electron 
collision frequency. The resulting two-dimensional equations express the 
conservation of electron flow and zero divergence of the total current. 

In a purely two-dimensional model everything is uniform parallel to the 
magnetic field and there is no flow of either species in the vertical direction. 
However, the same two-dimensional equations also represent the limit of 
vanishingly small parallel resistivity in which the electrons are allowed 
to flow freely up and down the magnetic field lines but ions are restrained 
to flow perpendicular to the magnetic field in response to the transverse 
electric fields. 

Figure II. 1 shows a schematic of the cross section of the ion cloud 
in the plane perpendicular to the magnetic field. The ion cloud is moving 
to the right with respect to the neutrals. In the two-dimensional model, 

C / SL 

the ratio of height-integrated conductivities; 2^/2^ , can be represented 
as the ratio of the cloud center density, N c> to the ambient ionospheric 
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Fig. II. I Simple Two-Dimensional Model 
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density, N q . Due to the slow rate of development of the ion cloud and the 
small plasma pressure relative to the magnetic field pressure (low (3), the 
ambient electric fields are electrostatic as expressed by the condition 

V x E = 0, Since the electrons do not experience any collisions they flow 

^ 2 

at the local IT x B/B velocity. The flow of the ionization is thus incom- 
pressible due to the assumption of electrostatic fields. The fact that the 
flow transverse to the magnetic field is incompressible is extremely 
important in theoretically determining the behavior of barium ion clouds. 
This assumption precludes the possibility of an electron density increase. 
The ion motion closely follows the electron motion but the ions do have a 
component of velocity in the direction of the electric field as a result of 
collisions with the neutrals. At high altitudes the ion velocity with respect 
to the neutrals is given approximately by 



?x B 1 IT 

B 2 *i B ' 


(II. 1 ) 


where /c. = m./eB T and T is the ion- neutral collision frequency. 

The simple set of coupled nonlinear differential equations that 
describe the motion of barium ion clouds transverse to the magnetic field 
are: 


9 N E x B 
9t + b 2 


VN 


0 ; 


(II. 2) 


V • Jp oc V • NET = 0; 


(II. 3) 


v x IT = o. 


(II. 4) 
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Equation (II. 2) is the continuity equation for the incompressible flow of the 

electrons at the local IT x B/B velocity. Equation (II. 3) describes the 

conservation of the current flow perpendicular to the magnetic field in 

the ionosphere and Eq. (II. 4) is the quasistatic assumption of the low- (3 

approximation. The Pedersen current, J , is proportional to the product 

P 

of the conductivity, N (proportional to the ion density), and the electric field, 

eT, all measured in the frame of the ambient neutrals. 

These equations have been used to investigate the stability of the 

3 4 5 

backside of barium ion clouds ’ , as well as other edges , and the non- 

linear development of striations^. Here we shall discuss how these simple 

equations have successfully been applied to describe the motion of barium 

7 6 

ion clouds , their deformations including the steepening of the backside , 

g 

and the predominant scale length of striations . First we describe in a 
simple way the relationship of this two-dimensional model to a more realistic 
layered ionosphere. 

When the barium ion cloud is considered to be of finite extent 
parallel to B , currents flowing in the ionosphere below the cloud are 
affected by the presence of the ion cloud and give rise to currents flowing 
parallel to the magnetic field. In Fig. II. 2 we indicate qualitatively the 
effect that the presence of the highly conducting barium ion cloud has on 
the current flow pattern in the ionosphere. Figure II. 2a is a schematic 
diagram in the plane perpendicular to the magnetic field showing the 
distortion of the current flow due to the presence of the more conducting 
ion cloud. In a completely two-dimensional model this same current flow 
pattern would exist at every altitude parallel to the magnetic field and drift 
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Fig. II. 2 Changes in the Ionospheric Current Flow Due to the Presence 
of a Highly Conducting Barium Ion Cloud; a) Height-Integrated 
Current Density Flowing Perpendicular to the Magnetic Field; 
This Same Pattern Results in the Cases of b) a One -Layered 
Ionosphere, c) a Two-Layered Ionosphere, and d) Local Plus 
Conjugate Ionosphere 
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as a "current cell" with the IT x B/B velocity. That situation is shown 
schematically in Fig. II. 2b which represents a cut in the yz plane through 
the center of the ion cloud. Note that the current density in the ion cloud 
is larger than in the less dense ionosphere. From Fig. II. 2a we see that 
this increase results from the concentration of the current in the ion cloud. 

In Fig. II. 2c the horizontal currents are the ion Pedersen currents 
whereas the parallel currents are due to electrons flowing along the magnetic 
field. Figure II. 2d represents a generalization to a three -layered iono- 
spheric model. The more conducting layer at the top of this figure 
represents a conjugate ionosphere which may be highly conducting and 
therefore can contribute a large perpendicular current density. Again the 
presence of a conducting ion cloud in the local ionosphere can effect the 
current flow in the conjugate ionosphere since electron currents flow freely 
along the magnetic field through the magnetosphere. 

While it is recognized that the diagrams given in Fig. II. 2c and II. 2d 
are qualitatively correct insofar as they indicate current flow in different 
layers of the ionosphere, they are strictly correct only if (1) there is no 
distinction between the barium ions and the ambient ionospheric ions; 

(2) if the neutral density is the same in all three layers; and (3) if the electric 
field which drives the currents is externally imposed on the complete system. 
The actual interaction of several different layers in the ionosphere with a 
barium ion cloud is only now beginning to be considered and is a very 
complex problem. It is clear that the plasma density will be increased 
where the electrons flow into the lower ionosphere and decreased where 
electrons leave to neutralize the charge accumulation in the cloud. Thus 
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the ionosphere subtended by the cloud will be locally perturbed and this 
perturbation will move with the cloud. 

We can see qualitatively that the effect of an ion cloud on the current 
flow pattern depends sensitively on the cloud conductivity compared to the 
conductivity of the ionosphere. The appropriate parameter that character- 
izes the cloud-ionosphere interaction is the ratio of the height-integrated 
Pedersen conductivity of the barium ion cloud to the height-integrated 
Pedersen conductivity of the ambient ionosphere. 

B. MOTION OF BARIUM ION CLOUDS 


Equation (II. 3) indicates that at any instant in time the electric field 
depends on the distribution of ionization in the ion cloud. In particular the 
polarization electric field due to the presence of the ion cloud reduces the 
total field at the center of the cloud as indicated at the top of Fig. II. 3. 

Asa result, the center of the ion cloud moves with a slower velocity with 
respect to the neutrals than does a background ion that is flowing with the 
velocity U q . The value of the electric field at the center of the ion cloud 
depends on the ratio of the height -integrated Pedersen conductivity of the 

Q 

ion cloud, 2 , to the height-integrated Pedersen conductivity of the back- 

cL 

ground ionosphere, 2^, as well as on the distribution of ionization in the 

ion cloud. For a uniform ionosphere and a constant neutral density, 

2p/2p can also be expressed as the ratio of height-integrated cloud density to 

ambient density, N /N . The electric field at the center of the ion cloud, 

7 c o 

and thus its motion, can be expressed in terms of a coupling parameter, t , , 
defined as the ratio of the electric field at the center of the ion cloud (in the 
frame of the ambient neutrals) to the external electric field, 
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Fig. II. 3 Effect of Polarization Field Induced by the Ion Cloud 


- 10 - 



(II. 5) 


? 


? = ?, 


pc 



where the subscript c represents the value at the center of the ion cloud. 
Since the difference in velocity between the ion cloud and the neutral wind 
is frame independent and equal to 1?^ x B/B , Eq. (IL 5) implies 




(II. 6) 


This equation also defines the frame invariant velocity of the ion cloud, 
U . It follows that £, ranges between one and zero for very small and 
very large ion clouds respectively; 


£ - 1 , 


V. 

x 


1? x B 
a 

B 2 


(N << N ); 
' c o" 


l - o. 



(N > > N ). 
c o 


(II. 7) 


Very small clouds or low density perturbations provide a direct measure 

of the ambient electric field, lT a = - x B , while very large clouds are 

completely tied to the neutral wind. 

Theoretical models have been explored that give analytic expressions 

for the dependence of t, on N c /N q . An example will be given in the next 

section. A well-known example corresponds to a cloud of constant density 

N with an ellipsoidal cross section as illustrated in Fig. II. 3. Two 
c 

interesting limits, a circular cloud with a = b and an elongated cloud with 
b << a are also shown. The former is the shape used byHaerendel, Lust, 
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and Reiger i n deriving their now classic equation for determining the 
ambient electric field from the motion of barium ion clouds. In the high 
altitude limit (k^ >> 1), their equation can be brought into the form of 
Eq. (II. 6) if the identification t, = 2/(1 + A') is made where X = (N q + N )/N . 
The elongated shape is most tied to the neutral wind and its motion corresponds 
to that of a field-aligned sheet of ionization. The degree of coupling also 
depends on the density profile of the ionization for large density clouds. A 

2 " 4/3 

cylindrically symmetric cloud with a density profile N oc (1 + r ) is 
less coupled to the neutral wind for the same value of peak column con- 
centration as indicated on the right-hand side of Fig. II. 3 and discussed in 
the next section. 

Figure II. 4 is a plot of the coupling parameter as a function of the 
ratio of cloud conductivity to ionospheric conductivity for several barium 
releases sponsored by ARPA in Project Secede. The estimated uncertainties 
in the values for £, are indicated by the vertical error bars, the horizontal 
error bars indicate a factor of two uncertainty in the estimates for the ratio 
of integrated conductivities. The two curves illustrate the two extremes 
in coupling that depend on the configuration of the ionization. The sheet- 
like shape is strongly coupled to the neutral wind, whereas the least coupled 
is the cylindrically symmetric cloud with a monotonic, radial density profile. 

A first glance at the data would suggest that the simple theory is inadequate 
to describe the degree of coupling of ion clouds to the neutral wind. We see 
that only the daytime and stammer morning releases come near to agreeing 
with the theory. Note the large disagreement for the winter evening releases 
since the scale on the bottom is logarithimic and covers several decades. 
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0.1 1.0 10 100 1000 
C 8767 Nc/No 

Fig. II. 4 Coupling Neglecting the Influence of the Conjugate Ionosphere. 

The open circle is event Quince while the half-filled circles 
are events Apple and Dogwood. In order from largest coupling 
parameter to smallest, the solid dots are the Plum Puff, 
Redwood, Plum, Spruce, and Nutmeg clouds respectively. 
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Thus the ion clouds are observed to separate from the corresponding neutral 
clouds much faster than the simple theory would allow. They behave as 
though the background ionosphere were far more conducting than the 
measurements of the ambient ionosphere indicate. 

The observation that the winter evening releases are those that are 
in most serious disagreement with the theory led to the suggestion that the 
conjugate ionosphere in the southern hemisphere was contributing to the 
background ionospheric conductivity. This suggestion seemed plausible 
because the ionosphere conjugate to Eglin (at 54°S, 103°W) is further 
west and was experiencing summer and hence was sunlit. On the other 
hand, the ionosphere conjugate to the summer morning releases (at 49°S, 
66°W) was experiencing winter and thus was still in darkness. More 
quantitatively, when the sun was 9° below the horizon after sunset at Eglin 
in the last half of January, 1971, and 9° below the horizon before sunrise 
at Puerto Rico in the first half of May, 1969, the sun was 26° above the 
horizon and 2 3° below the horizon, respectively, at the corresponding 
conjugate ionospheres. Hence it is reasonable to assume that the height- 
integrated Pedersen conductivity of the conjugate ionosphere for the winter 
evening releases had a typical daytime value of 5. 0 mho. By adding these 
conductivities to the local values, we obtain new estimates for the ratio of 
cloud to ionospheric conductivity. 

Figure IL 5 is a replotting of the data assuming that the conjugate 
ionosphere does contribute to the background ionospheric conductivity. By 
making this assumption we see that there is no longer any disagreement 
between the simple theory and the data. It is clear that these clouds do 
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Fig. II. 5 Coupling Including the Influence of the Conjugate Ionosphere. 

The dots represent the same clouds shown in the preceding 
figure. The solid square represents the point obtained from 
numerical modeling. 
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not move as though they were shaped like sheets. However, the data do not 
clearly distinguish between the two other cylindrical models whose theoretical 
curves are shown. Zabusky, Perkins, and Doles ^ have solved Eqs. (H. 2) 
through (II. 4) by finite difference techniques for a cloud with a cylindrically 
symmetric initial density distribution, N = N q [ 1 +4(1 - r^/4)^] for r < 2. 
The velocity of the peak density in the E x B frame is 0. 63 V . The 
corresponding value for the coupling parameter, £, , is 0. 37. This value 
for N^/N q = 4 is shown as a solid square in Fig. II. 5. We note that this 
point is in excellent agreement with the analytical models, lying closer to 
the curve with a density profile (which should provide an upper limit) than 
to the curve for a cylinder with constant density. 

We note that the agreement with the theoretical model includes data 
in which the ratio of conductivities ranges over two and one -half orders of 
magnitude. This agreement is highly suggestive that the conjugate ionosphere 
must participate in the dynamics of the motion of barium ion clouds. We 
estimate that an Alfven wave (the propagation velocity of an electromagnetic 
disturbance) can travel from the northern to the southern hemisphere 
and back in half a minute, so there is sufficient time to communicate the 
polarization electric fields between the two hemispheres. The most 
sensitive test of the theory would be to make a 12 kg size release at dawn 
and dusk. The degree of coupling with the neutral wind should change from 
about 65% (£, = 0. 35) to about 30% (C, = 0. 7) for the assumed change of a 
factor of five in the ratio of conductivities. 

There are two other possible explanations for the lack of agreement 
between data and theory shown in Fig. II. 4; either the data are plotted 
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incorrectly, or the theoretical treatment is wrong. The largest uncertainty 
in the velocities listed is the ambient ionospheric velocity, IT x §/B , which 
has been assumed to be less than 20 m/s where better estimates were not 
available. The remaining possibility is that the theory does not correctly 
describe the coupling of the barium cloud to the rest of the ionosphere. It 
is true that the complications of this interaction are hidden by ascribing 
to the ionosphere a height-integrated conductivity, a procedure that would 
be strictly valid only if the value of outside the cloud were constant and 
the same as for barium. There is no assessment yet as to the effect to be 
expected from a more rigorous treatment of the dynamics of the rest of the 
ionosphere. On the other hand, the assumption that the reconciliation 
between data and theory is achieved by including the effect of the conjugate 
ionosphere is a simple and plausible one. 

There are several important implications that follow from assuming 
that the conjugate ionosphere plays a role in the ion cloud dynamics. First 
is that the cloud -produced polarization electric fields must project into the 
southern hemisphere which suggests that the assumption that the earth's 
magnetic field lines are equipotentials for a large distance is reasonable. 

A second conclusion is that there must be magnetic -field -aligned currents 
in the magnetosphere, and a third conclusion is that the presence of the ion 
cloud and its resulting motion in the northern hemisphere must cause a 
disturbance in the ion density in the conjugate ionosphere in the southern 
hemisphere. An independent verification of any of these effects would be 
very valuable. 

Let us adopt a simplified model of an ion cloud in the ionospheric - 
magnetospheric system as shown in Fig. II. 6 in order to make a quantitative 
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estimate of some of these effects. We find that the magnetic field lines are 
highly conducting so that the electric fields do not suffer any appreciable 
resistive losses in being projected into the conjugate ionosphere. The 
most important effect is the disturbance in electron density that the current 
flow in the background ionosphere is likely to make. The parallel currents 
must be closed by Pedersen currents which are due to the flow of the 
background ions. This ion flow will lead to significant changes in the local 
ionospheric density. These changes should be observable by ionosondes. 

In particular, for the winter evening releases, the motion of the ion clouds 
should have produced large changes in electron density in the conjugate 
ionosphere in the southern hemisphere. These changes should have been 
detectable had ionosondes been deployed there. 

C. DEFORMATIONS OF ION CLOUDS 

The electric field varies over the distribution of density in the ion 

cloud as discussed above and as shown schematically in Fig. II. 3. Since 

, , o 

the ionization moves with the local E x B/B velocity, different parts of the 

ion cloud move with different velocities. Thus, as the ion cloud moves, its 

shape distorts. Here we shall present the results of an analytic investigation 

of the changes in shape of equal density contours obtained by expanding the 

solution of the nonlinear equations for the equipotentials (streamlines) as a 

power series in time. We shall see that there are qualitative differences 

between the contours obtained here and thos recently presented by Simon 
17 

and Sleeper. We shall point out that they obtained incorrect deformations 
because they did not carry through their expansion procedure correctly. 


We write the electric field as the sum of the externally imposed 
uniform field, £" q , plus the polarization field, IT^, due to the presence of the 
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ion cloud; 


r = r o + r p = F o - * * p ( u -8) 

Equation (II, 3) is conveniently rewritten 

^ 2 <^p = {^" 0 " V <f>p) ‘ V Li (II. 9) 

where 

L = in N/N q = in ( 1 + j^-\ . 

AN is the additional density due to the presence of the ion cloud. 

Equation (II. 9) determines the spatial variation of the polarization field 
as a function of the ion cloud density distribution. Once cf> is known, 

Eq. (II. 2) can be used to obtain the changes in density that result. The 
new density distribution including these changes then produces a new 
polarization field distribution. It is this last step that was omitted by 
Simon and Sleeper^ resulting in their erroneous deformation contours. 

We have found a useful self-consistent solution of Eq. (II. 9): 


po 


_ b r cos 0 

-- P E 0 

b + r 


(II. 10) 


Lq(z) = in 


1 + 


AN 

N 


= in 


(1 + z) 


z +(2 + p)z + (l-p) 


Va 


in 


2 z + 2 -f p + 


V8p 


±JL 


ip + p 


2 z + 2 + 


p - V8p 


+ P 


(II. 1 1 ) 


where z = ( r/b ) . According to this solution, the polarization electric field 
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at the center of the cloud is 


\ = -P^ 0 - (H.12) 

Hence, the parameter p represents the amount of reduction of the ha ck- 
ground field at the center of the cloud due to the cloud’s presence and is 
equal to 1 - t, where £, was defined in the previous section. This decrease 
in the electric field at the cloud's center makes it easy to see that the region 
of maximum density does not keep up with the background ionization and 
thus falls back toward the neutral cloud causing the backside to steepen. 

From Eq. (II. 11) , we can determine the cloud center density as a 
function of the parameter p by setting z = 0: 


L 

o 


( 0 ) 


in 



in 


1 -p 


Vs, 


in 


2 + 


V So 


p + p 


4(1 - p) 


J 


(II. 13) 


By expanding this equation in the limits p << 1 and p — 1, we obtain the 

relation between p and the cloud-center density for low- and high-density 
clouds ; 


P 


( N /2N 
I c 7 o 



4E 


<V N o> 


374 


N « N 
c o 


N » N 
c o 


(II. 14) 


The half-width, a, at half-max of the cloud is determined by the equation 
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In the limit of small p (low cloud density), Eq. (II. 11) can be 
expanded to yield 

AN 2 p 


N 


2 2 

(d + r ) 


r ; (N c « N o ) 


for which the half-width, a, is given as 

a = (yjT - l) l / Z b = 0.644b. 

In the opposite limit of large cloud density, we may obtain 


(II. 16) 


(II. 17) 


AN 

N 


2/3 
2 4 /3 


; (N » N ) 
c o 


(1 - p + 3r ) 


(II. 18) 


for which the half-width, a, is determined as 


a ~ 


(2 


3 / 4 . n 1^-fi 


-,1/2 


b - 


0 . 626 


© 


3/8 


(II. 19) 


The normalized cloud density profiles AN(r)/N c> as determined by Eq. (II. 11) 
for the two limiting values of the parameter N c /N q and normalized by their 
half-width, a, are shown in Fig. IL 7. A normalized gaussian profile is also 
shown for comparison. It is seen that the shape of the cloud is fairly 
reasonable, falling off at large distances as r 4 and r for the two 


cases shown. 




Fig. II. 7 Normalized Density Profiles in the Limits of Very Low 
Density and Very High Density Clouds Compared to a 
Gaussian Profile 
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In the frame moving with the velocity of the background ionosphere, 
U q = x B/B , the equation that describes the deformations of the density 


contours can be written 


9 N 

at 



x B 
2 


VN 


(II. 20) 


From this equation it is easy to see that the equipotentials determined by 
= constant are streamlines of the flow. An example of the streamlines 
derived from Eq. (II. 10) are shown in Fig. II. 8. From Eq. (II. 20) we may 
calculate the form that the initial displacements within the cloud take. Since 
the cloud is initially cylindrically symmetric, we need only consider the 
radial components of the displacements of the contours. We thus have 
for a small, initial time step 


i w 'r 

A r = v At = — =j- a - At = g(r) sin 6 At (II. 2 1 ) 

r ro a 

where g is a function of radius alone. Hence each equal density contour, 
originally a set of concentric circles as shown in Fig. II. 8, is displaced 
initially to the left. There is no change in shape of each density contour. 

If the streamlines shown in Fig. II. 8 were constant in time, then 
the changes in shape of the equal density contour could be easily determined 
by following the streamlines shown. The deformations that would result 
are shown in Fig. IL 9a. In actuality, the streamlines are not constant in time 
but are governed by the time -dependent Eqs. (II. 9) and (IL 20). The changes in 
density resulting from the displacements given by Eq. (11.21) produces a 


- 24 - 




C 7568 


Fig. II. 8 


Equipotentials Resulting From the Solution Given in 
Equation (II. 10) 







new contribution to that changes the streamlines. We shall calculate 
these changes to next order in time in the form 


(r, 0, t) = 



(r, e 


) 


1 + pt fj (r, e 


) 


(II. 22) 


L (r, 9 , t) = L q (r) + pt Lj (r, 9, p) 


(II. 23) 


where <j> and L q are given by Eqs. (II. 11) and (II. 10), respectively. L.^ 
is determined by Eq. (II. 20) and the functions and L q . can then be 

used in Eq. (II. 9) to determine f . It is not possible to solve the resulting 
equation for f^ exactly for arbitrary p. The resulting equation is expanded 
in a power series in p and the first terms of that series have been determined. 
The resulting expressions are: 


Lj (r, 0 , p ) = — 2 2 

b + r 


b- 9 p ] 


9 r 


sin 9 ; 


(II. 24) 


f 2 r . „ 

f l = - 3 b Sln 9 


1+ ^ 3 + Z z - IT ® — 4 + 0(p 2 ) ] (IL 2 5) 

(1 + z) (1 + z) J 


2 

where z = (r/b) . It is found that the term in f^ which depends on p makes 
negligible difference in the resulting deformations of density contours even 
for the value p = 1. Hence, the neglect of this term altogether is probably 
reasonable except for the very highest density clouds. 

We now have a solution for the time dependent stream functions. The 
resulting deformations in the density contours when they are allowed to follow 
these time dependent stream functions are shown in Fig. II. 9b. Note that 
the higher order term, f ^ , makes a qualitative difference from the changes 
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in shape shown in Fig. II. 9a. Note that the contours are convex outward 
on the righthand or front side of the barium ion cloud when the more 
accurate expression is used. We also note that the backside has become 
steeper during the same time interval. 

The deformations in Fig. IL 9b show several interesting features 
for which we now have an analytic handle. The first is the steepening 
of the backside which results from the polarization electric field, IT , 
due to the density gradients (or, alternatively, to the increased drag at 
the center of the ion cloud which pulls the center in the direction of the 
neutral cloud). The second is the increase in width in the direction 
transverse to its motion as shown by the high-density contours near the 
center of the cloud. Since the flow described by Eq. (II. 20) is incompressible, 
the area within each equal density contour is constant. This increase in 
width observationally would appear as though the cloud were diffusing and 
thus this model which assumes no diffusion provides an explanation 
for the anomalously large "diffusion coefficient" which has been reported 
for diffusion of ion clouds perpendicular to B. We can also see that there 
is an elongation of the lower density contours in the direction of cloud 
motion. In particular, at times longer than that for which Eq. (11.22) is 
valid, a long thin neck forms where the center of the cloud had been. 

Lastly, the deformations shown provide an explanation of the fact that the 
trailing edge of ion clouds develop a large density gradient when seen from 
a variety of views. 

We shall not discuss the solution for the diffusion and motion of 

12 

barium ion clouds derived earlier by Simon. Both results differ from 
the model presented here and are valid only for the case N c « N q . We 
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restrict our attention to the high altitude situation » 1) where these 

qualitative differences are quantitatively small. We note that in this limit 
of negligible diffusion Simon and Sleeper's^ Eqs. (8) and (16) are 
equivalent to our Eqs. (II. 2) and (II. 3) for the dynamics of barium ion 
clouds. However, they have neglected the contributions to the deformations 
resulting from the changes in density. Thus their deformed contours (see 
their Fig. 1) correspond to the incorrect one shown in Fig. II. 9a. The 
reason for their error is probably related to the fact that they attempted to 
present a solution correct to first order in p. The term that gives rise to 
the differences is the second term of Eq. (II. 22) which is of second order 
in p for any finite time. However, for times that are appreciable fractions 
of Ba/pE Q , this term gives rise to the qualitative differences shown. 

The deformations given by Fig. IL 9b agree qualitatively with 
observations of barium ion clouds while those represented by Fig. II. 9a 
and given by Simon and Sleeper do not. The most striking feature is the 
steepening of the backside (the bunching of the density contours on the left 
hand side). The broadening of the highest density contours transverse to 
the direction of motion is also evident as well as the tendency for the lower 
density contours to elongate in the direction of motion. Zabusky, Perkins, 
and Doles^ have numerically integrated Eqs. (II. 2) and (II. 3) for the case 
of an initially cylindrically symmetric ion cloud with N /N q = 4. Their 
results agree qualitatively with that shown in Fig. II. 9b. 

D. SCALE SIZE OF STRIATIONS 

3 

The major criticism of Linson and Workman's early work is that 
it represents the stability of an infinite half -plane with infinite density far 


- 29 - 



to the right and zero density far to the left. That analysis has recently been 
extended in several ways, primarily by taking into account the finite extent 
of the back edge of an ion cloud as well as the finite ratio of ion cloud 
conductivity to the background conductivity. 

The basic equations used in the stability analysis are Eqs. (II. 2) and 
(IL 3) with two modifications. The first is that if we include electron 

collisions, we obtain a term V- D | V N on the right-hand side of Eq. (II. 2). 

2 2 

Here D x = K(T g + T.)m e v J e B is the ambipolar electron-collision-controlled 
diffusion coefficient. The electron collision frequency, v g , is taken to be 
a constant even though for high density clouds collisions with ions dominate 
collisions with neutrals. The second modification is related to the fact that 
we will be analyzing a one -dimensional slab model as an approximation to a 
local analysis on the back edge of an ion cloud. Hence the unperturbed 
electric field that should appear in Eq. (II. 2) is the electric field E^ = 1, 1?^ 
that exists at the edge of the ion cloud due to its polarization. Alternatively, 
the unperturbed velocity E x B/B should be the velocity of the back edge of 
the ion cloud, tT c , with respect to the neutrals given by Eq. (II. 6). For large 
clouds, this velocity is considerably different from the ambient flow velocity 
U q = E q x B/B 2 (see Fig. II. 5). 

The zeroth order density profile whose stability is analyzed is 

x < 0 

0 < x < x = d £n (1 + N /N ) (II. 26) 

o do 

X < X 
o 


N 


N = { N e' 

o 


c/d 


N + N 
o c 
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Note that the logarithmic density gradient length d = x q /£ n (1 + a ) is given 
by the extent of the back edge , x , and the ratio of cloud to ionospheric 
conductivities, a = N /N q . A straight-forward linearization of Eqs. (II. 2) 
and (II. 3) results in the following equation for the x-variation of the perturbed 
potential, , 


it ii 

N (f> + (j) N 




0 


(II. 27) 


where the growth rate y is given by 


V 


k 2 D 


JL‘ 


(II. 28) 


The y-variation of the perturbations has been taken to be exp(i k y). The 

problem is solved when the boundary conditions of vanishing potential at 

x — - _+_ oo determine unique values of y 

Equation (II. 27) can be solved piecewise for the unperturbed density 

r 

distribution given by Eq. (II. 26). By applying the continuity of </> and <j> at 
x = 0 and x = x q , and defining 


K = 



(II. 29) 


the equation that determines y implicitly can be put in the form 


tank Kx 

o 


2 K 



(II. 30) 
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Some useful limiting solutions of this equation are: 


^ kU tanh -k 2 D, ; 

c 2 _L 


U 


2 / 2 

l _ 1 + it / a 

4(kd) 2 


k 2 D, 


V = 


a kU 


2 + a kd 


-k 2 D, 


k U 

2 2 l/ 2 

V. (1 + 4k d ) -kd 


-k ; 


kd << 1 


kd » 1 


(II. 31) 


a << 1 (d fixed) 


a » 1 


From these forms we see that long wavelengths are stabilized by the finite 
size of the cloud while short wavelengths are stabilized by diffusion. There 
is a wavelength for maximum growth that can easily be derived from the 
approximate expression 

U 

c 


1 - 


- k 2 D. 


4(kd) 


(II. 32 ) 


yielding 



(II. 33) 


A typical linear -wave spectrum is shown in Fig. II. 10. There is a 
wavelength for maximum growth, but it should be pointed out that the 
wavenumber scale is logarithmic and there is a significant growth rate for 
a wide range of different scale sizes. It is reasonable to expect a variety 
of wavelengths attempting to grow initially. It should be emphasized that 
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Fig. II. 10 


Linear Wave Spectrum 





the value of the wavelength for maximum growth depends sensitively on the 
parameter d = J V £n N | * which so far can only be estimated for large density 
barium ion clouds that have steepened. 

We will now compare the results of the above analysis with observa- 
tions of barium ion clouds made at Eglin AFB in January 1971. The unique 
aspect of these releases is that three clouds, Nutmeg, Spruce, and Redwood, 
released at altitudes of 150 km, 185 km, and 255 km respectively, passed 
through the magnetic zenith at a time when they were striated. Figs. II. 11 
and II. 12 are photographs of Redwood and Spruce, respectively. We note 
that there is a good deal of structure evident as the striations are highly 
nonlinear. However, there is an overall sheet-like structure evident in 
both photographs. The width of the sheets are compared to the results of 
the linear theory in Table II. 1. 

TABLE II. 1 Dominant Scale Size of Striations 



The appropriate values for U c and were used for each ion cloud. We see 
that the uniform assumption that d = 500 m does not give good agreement 
with the observed dominant wavelength. However, if d = 500 m is assumed 
appropriate for Spruce and d is scaled to the initial transverse width, a, of 
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Fig. II. 11 Redwood Striations Viewed Driectly Up the Magnetic Field 
Line (Technology International Corporation) 
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Fig. II. 12 Spruce Striations Viewed Directly Up the Magnetic Field 
Line (Technology International Corporation) 
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each of the other two clouds, then the wavelength for peak growth gives 
excellent agreement with the observed dominant widths. 
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III. ION STREAMING INSTABILITIES WITH APPLICATION 
TO COLLISIONLESS SHOCK WAVE STRUCTURE 


A. INTRODUCTION 

The instabilities that result from the presence of two interpenetrating 
ion streams immersed in an electron background may provide the mechanism 
for collisionless momentum coupling between the two streams. The presence 
of a weak magnetic field can have a significant effect on the physical process. 
First, the magnetic field can inhibit the electrons from shorting out ion 
plasma oscillations for wavelengths long compared to the electron gyro- 
radius for propagation perpendicular to the magnetic field; second, gradients 
in the magnetic field give rise to electron currents that can drive ion acoustic 
waves unstable and increase the effective electron collision frequency; and 
third, whistler mode -- ion beam interactions are likely to be important, 
and the existence of whistler modes depends on the presence of a magnetic 
field. 

The existence of collisionless plasma shock waves has provided the 

impetus for investigating mechanisms that can provide the anomalous 

dissipation that is required. Many experiments have shown that plasma 

turbulence is associated with these shock wave structures. Numerous 

linear plasma instabilities have been proposed as the mechanism that 

generates the observed turbulence. An overall review of the status of 

theoretical and experimental work in collisionless shock waves has been 

given by Friedman, et al. ^ We shall briefly review here the status of 

1 /2 

perpendicular, (ir/ 2) - (m e /rm) ' < 0 < 7r/2, (© being the angle between 
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the flow direction and magnetic field) and oblique, 0 < © < [tt/Z) - (m e/ m i } 

shock wave structures together with the different streaming instabilities and 

their possible roles in the formation of these structures. We shall then 

present an analysis of a new configuration leading to instability with possible 

application to parallel, intermediate Alfven Mach number shocks. 

For low (3, perpendicular shocks, when > 5-10 T., a relative 

drift between the ions and electrons in excess of the ion acoustic speed can 

2 3 

excite unstable ion acoustic waves. ' The ensuing turbulent fluctuation 
2 

spectrum <l€ k l > stochastically heats the electrons and the shock 
structure is controlled entirely by the electron resistivity associated with 
< | E^ | >. This kind of shock structure can be oscillatory or monotonic 

(depending on the amount of turbulent dissipation) and it exists as long as the 
downstream plasma flow velocity, V^, is greater than or equal to the down- 
stream ion acoustic speed, Cg^. When < Cg^ (corresponding to an 

upstream Alfven Mach number M^ ^ > M = 2. 76), a narrow gasdynamic 

4 

subshock is embedded in the magnetic shock layer. At these higher 

upstream flow velocities, it is the presence of ion viscous dissipation 

dominating (only) the subshock structure that provides the requisite 

degradation of the additional directed energy into thermal energy. 

5 

Recent studies by Papadopoulos et al suggest that two ion beams 
counter streaming in a background rest frame of cool magnetized electrons 
can also effectively couple momentum from one stream to the other at 
characteristic frequencies and growth rates comparable with the lower 
hybrid frequency. This instability is operative only if the Alfven Mach 
number, M^, for the relative drift between the ion streams satisfies the 
inequality 
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(III-l) 


M 


< 2(1 + P e ) 


1/2 


and if the streaming direction is nearly perpendicular to the magnetic field. 

We shall introduce the concept that the relative drift is the difference between 

the upstream and downstream velocities, and V^, respectively, i. e. = 

Vj - V^. If we define the Alfven Mach number in terms of the local value 

1 /2 

of the Alfven velocity which varies as B/n ' , we can express the limitation 

implied by Eq. (Ill - 1 ) in terms of the upstream Mach number, M^ = Vj/C^. 
For instance, in the case of a low (3, perpendicular shock structure with 
the interaction region taken to be the center of the shock layer so that 
(3(x) ~ 1, n(x) ~ (l/2 )(n^ + n^), it can be shown that 

1/2 


2 (2I 1 / 2 > M d (*) = M A1 


(III-2 ) 


where j - n^/n^ is the density compression ratio. By assuming y = 5/3, 

T may be related to the shock Alfven Mach number M^ through the relation 


^25 + 52 M a } + 4M a , 4 - I - M 


A 1 


A A1 


A1 


(III- 31 


valid only when the unshocked plasma is cold. Thus, from Eqs. (III-2) and 

(III- 3) this magnetized ion-ion instability can operate only for M^ < 6. 

6 7 

The "modified" two stream instability ’ has also been proposed for 

nearly perpendicular shocks. This instability is driven by relative streaming 

of ions and electrons across a magnetic field. It is so named because its 

dispersion relation, in the fluid limit, is similar to that for the Buneman 

8 

two-stream instability . Unlike the ion acoustic instability which selectively 
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heats electrons, the "modified" two-stream instability produces large 
and comparable ion and electron heating. Like the magnetic ion-ion 
interaction discussed in the previous paragraph, the characteristic frequency 
and growth rate for the modified instability is comparable to the lower 
hybrid frequency and it is operative only if the Alfven Mach number for the 
relative drift between the ions and electrons satisfies the inequality, 

M d < < ! + Pe ,l/2 ' 

Significant among recent advances in the area of magnetic streaming 
instabilities and their applications to oblique shock structures in collisionless 
plasmas are the beam -whistler (BW)^ and beam-cyclotron (BC)^’ ^ 
instabilities. These are characterized by strong interactions between one 
of the two counterstreaming ion acoustic beam modes and the whistler or 
electron cyclotron mode (see Fig. Ill - 1 ). While the BC instability does 
result in electron heating, it nevertheless is inadequate as a mechanism for 
ion-ion momentum coupling. Here the electron cyclotron modes are heavily 
damped and do not interact with the beam modes. On the other hand, BW 

i /z 

instabilities, operative in the range 0 < © < ( 7r/ 2 ) and < (rm/m^) ' X 

cos © , do result in turbulence which can suitably couple the two ion beams. 
Clearly, the application of BW interactions to oblique shock structure is 
expecially attractive since they are operative over a substantially larger 
range of shock Mach numbers than are the magnetic counterstreaming ion- 
ion and modified two stream instabilities discussed earlier. 


We now briefly review the oblique collisionless shock structure study 

q 

of Lindeman and Drummond. Their analysis points to the existance of four 


critical Mach numbers M , M , M 


The first of these, cited 


earlier for perpendicular shocks, is the value of the upstream Alfven Mach 


- 42 - 



U) 



Fig. IIL 1 Schematic of the Dispersion Relation Showing the Interaction 
Regions Leading to Beam-Cyclotron (BC) and Beam- Whistler 
(BW) Instabilities. 
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number for which the downstream flow velocity equals the ion acoustic 

speed. When ^ < M , the shock structure is a large amplitude 

whistler (in the rest frame of the shock front). The damping of the whistler 
and subsequent electron heating are brought about by the anomalous resistivity. 


When > M, a nonmagnetic subshock appears downstream of the whistler 

region. It was found, however, that this subshock is a single large amplitude 

ion acoustic wave of scale length of order of the Debeye length; this refutes 

the usual notion of the anomalous ion viscosity as a broadening mechanism 

for subshocks. When > M , a fraction of the incoming ions is reflected 

from the potential barrier of the ion acoustic wave. These reflected ions 

then interact with the counterstreaming upstream ions through a BW 

instability resulting in pickup and heating of the reflected ions by the incoming 

ions. Since the relative drift velocity between the incoming and reflected 

g 

particles s 2V, Lindeman and Drummond predict a third critical Mach 


number M 


si- si- sl- 

•¥• "i> t» 


j /o si- si- st- 

0. 5 (m./2m ) ' cos 0 , such that for M. . > M , the 


i' e' ’ - ' A1 

above BW instability which couples the counter streaming ions is inoperative. 

*** * . \/2 

Finally, when M > M =0. 5(rm/m e ) cos 0 , the whistler can no longer 

stand in the shock layer (in the rest frame of the shock) so that its structure 

is completely dominated by the ion acoustic mode. 

It is interesting to note that for parallel propagation (k (f ), the 

ion acoustic modes of the counterstreaming ions uncouple from the whistler 

mode so that in the assumed frequency range 0 ^ « | <*j | , | co - kU | « 

9 

fi of Lindeman and Drummond, the BW instability no longer exists. On 
the other hand, if one relaxes the frequency restriction to lower frequencies, 


i. e. 0 ~ 


co 


CO - kU « SI , then even in the case of parallel propa- 


gation a beam mode can, in fact, enter into unstable interaction with the 
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with the whistler leading to coupling of the two ion beams. This problem 
and its application to parallel shock wave structure are the topics of 
central interest in this report. While our analysis is restrictive inasmuch 
as we concentrate only on the k j | , V || B q case, it is nevertheless useful 

to know if ion dissipation can be generated. 

In Section ILL B, we analyze the linear instability of two counter - 
streaming ion beams in the presence of magnetized electrons; here, the 
analysis is carried out in the rest frame of the electrons. In Section III. C 
we analyze the resulting dispersion relation in the rest frame of the parallel 
shock at different upstream Mach numbers and at different stations in the 
shock layer. The result of this analysis is a map showing the range of 
Mach numbers as a function of position into the shock for which unstable 
waves can stand inside the shock. 

B. LINEAR THEORY 

We analyze here the instability arising in a collisionless plasma 
configuration consisting of two ion beams counterstreaming along the 
magnetic field (V g ||B o ) through an electron background. V g represents 
the streaming velocity of s LX1 species of ions. The ion and electron 
motions are modeled by fluid equations valid for wavelengths long compared 
with the electron gyroradius. Restricting ourselves to the case of parallel 
propagation (k || B^ ), the appropriate electromagnetic dispersion 
relation then has the form: 

,2 2 

V” " 6 11 w ^ = — ie 12 (III-4) 

00 

where, in the electron rest frame, 
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and 


e 12 (k, co 


^ i (0> 2 r , a /w 2 )(w - k • VJ J2„. 
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(III-5b) 


= Z | e | B^/rm, ^ ce = | e | -® 0 / m e deonoting the ion and electron cyclotron 

th. 

frequencies and COp S » CO the s ion beam (s = 1,2) and electron plasma 
frequencies, respectively. We note that within the fluid approximation, the 
dispersion relation Eq. ( III - 4 ) is unaffected by the presence of electron 
thermal pressure since the thermal velocity does not appear in the 
11, 22, and 12 components of the dielectric tensor in the fluid limit for 
the case of parallel propagation. 

We now make the low frequency approximations that 1) the displacement 

current is negligibly small (co « kc), and 2) the electrons have zero mass 

( CO I « , co ). Introducing the dimensionless densities 5 = Zn_/n , 

I ce pe 2 e 

(1 - 6) = Zn^/n^ (it is understood that the composite ion beams -electron 
plasma system is charge neutral), the dispersion relation simplifies to 


k 2 c 2 


CO 


(1 - 6 ) ( co - k 
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ci 


co - k • V + n 


Cl 


co - k 


v, + n . 

~2 — ci 
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(III-6a, b) 

where Eq. (III-6a) refers to the upper signs, and Eq. (III-6b) refers to the 
lower signs. We also apply the zero current constraint 
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0. 


(Ill - 7 ) 


u - 6) y 1 + 6 y 2 = 


For the low frequency case under consideration here, it is crucial that 
condition (III-7) be met for arbitrary beam strength (0 < 5 < l)(see Ref, 12). 
In the sequel, Eq. (III-6b) corresponding to the lower sign, can be dropped 
from further consideration since it is identical to Eq. (III-6a) with a J — - CO 
and k - k. It is convenient to introduce here the relative drift velocity 
U = V 2 - y r so that from Eq. (III-7), = -6 U, = (1 - 5)U, 

and Eq. (III-6a) can be simplified to 


2 2 
k c 


2 

Cl) 



ci 


(1 - 6) ft . 6ft 

ci ci _ Q 

(_t) + 6k » U + ft . a) - (1 - 6)k • U + ft 

~ ci — ~ ci 

(III- 8) 


Now, for 0 < 6 < 0. 5, it is clear that ion stream 1 has more the 
role of "plasma" while 2 has more the role of "beam". When 0. 5 < 6 < 1, 
these roles are interchanged. Since, however, the corresponding interchanges 


6 — * (1 - 6), k • U — -k • U 

leave Eq. (Ill - 8 ) unchanged, analysis of the range 0. 5 < 6 < 1.0 would be 
redundant. It is therefore sufficient to investigate only 0 < 5 < 0. 5. 

For given values of the parameter 5 and = U/C ^ (it is convenient 

to normalize U with respect to the Alfven speed ), our task consists in 

solving Eq. (III-8) for aj = a) (k); here k is considered to be real and c 0 
can be real or complex. Since Eq, (III- 8) is a cubic equation in ct> with real 
coefficients for real k, complex roots necessarily occur in conjugate pairs. 
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Hence, whenever a growing mode (Im c 0 > 0) is present, there is also a 
decaying mode (Im u> < 0). 

For 6=0, one recovers the usual whistler (W + ) and ion-cyclotron 
(IC + ) dispersion branches given analytically by 


CO 


= k C. /2 fl . + k C. (1 + 


ci 


k 2 C, 


740 J) 


1/2 


ci 


(III-9) 


For weak beams (6 « 1), the intersection points P, Q. R (see Fig. III-2) 
between the ion beam mode 


co = (l-fi)kU - a . 

ci 


(III- 10) 


and the W + and IC_ branches of Eq. (III-9) suggest possible strong interactions 

1 3 

near these points. Indeed Kovner's analysis shows that beam -IC_ inter- 
action is stable while the beam- W + interactions are unstable with growth 
rates 


. k U 
ci St . 


■i 2 ..2 _ 2 

k. U - fi 
a ci 


172 


(ft = P.Q). 


(Ill- 1 1 ) 


where, from Eqs. (III-9) and (III-10), kp and k^ depend on the slope (=M^) 
of the beam line. Note that the beam-cyclotron (BC) intersection shown in 
Fig. III.l does not appear in Fig. III. 2 due to our approximation m g 0. 

We have extended Kovner's weak beam study by solving Eq. (III-8)for 
arbitrary values of 6 (0< 6 <0. 5). By expanding the dispersion relation 
(III-8) for small k, one readily finds 

2 1 '/ 2 


Im w = kC 


Jg (1 - 6) M d 2 - 1 


(III- 12] 
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Fig. III. 2 Superposition of the Whistler (W + ) and Ion Cyclotron (IC + ) 
Modes for 6 = 0 with a Weak Bearm Mode. The intersection 
points R, P, and Q represent regions where strong interactions 
can occur. 
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so that these long wavelength modes are unstable only if 


1 

M d > [6 (1 - 6)] l / z 


(III- 13) 


This condition is thus sufficient for instability, but is not necessary as we 
shall see below. We should emphasize that due to the application of the 
condition (III-7), we are dealing only with zero-current instabilities and thus 
the stronger, current-driven instabilities often treated are excluded. It is 
easy to show that there must be complex roots of Eq. (Ill- 8 ) for values of 
k and U determined by setting u) given by Eq. (Ill- 10) equal to cu given 
by Eq. (III-9). The resulting relationship between U and k is 


£2 


U 


ci 


(1 - 5)k 


1 + V— + 


£2 


k 2 2 2 

(l+k Z C A 74fi c p 


ci 


£2 


. (Ill- 14) 


ci 


By appropriate manipulation, Eq. (III-8) can be brought into the form 


(co + 6kU + £2 . ) (co - u) )(u) - (jj , )(u> - (1 -6) k U + £2 . ) = -5 ( l - 6)S2 2 k 2 U 2 . 

(III- 15) 

where we have introduced a fourth, extraneous root, u) = - £2 . , in order 

ci 

to produce the above quartic equation in o> . When Eq. (Ill- 14) is satisfied, 

the third and fourth factors on the lefthand side of Eq. (Ill -15) are identical 

and hence positive for real u ) . The righthand side is always negative while 

the lefthand side is negative only for -5 k U - £2.^ < 0) < 0J_ < 0. In this 

range there is only one real root,, co, , such that -5kU - £2 . < co, < - £2 . , in 

addition to the extraneous root at - £2 .. Hence the other two roots of 

ci 
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Eq. (Ill- 15) and thus Eq. (III-8) are complex for k and U determined by 
Eq. (HI- 14). 


Dispersion curves obtained from Eq. ( III - 8 ) are displayed in 
Figs. III. 3a - III. 3d for values of 6 = 0. 1 and 0. 5 and drift Mach numbers 
= 3 and 7. Note that kU has been taken to be negative here as opposed 
to the positive value shown in Fig. III. 2. We find that for 6 = 0. 1, the peak 
growth rate increases only slightly from y ~ 0. 32 £2 at = 3 to 
y * 0. 38 at = 7, while for 6 - 0. 5, the peak growth rate increases 
more markedly with M^, rising from y = 0. 34 ^ at = 3 to y - 0. 72 SI ^ 

at = 7. Note that for 6 = 0. 5, the dispersion curves are symmetric in k 
as it is impossible to distinguish the beam from the plasma while for 
6 = 0. 1 the beam-like mode can be clearly identified. Also note that in 

Fig. III. 3a the region near k = 0 is stable in agreement with Eq. (Ill- 13) 
but that there is instability for a finite range of negative k. 

The regions of instability in k, U space for 6 = 0. 1 and 0. 5, obtained 
by setting the discriminant of Eq. (III-8) equal to zero, are shown in 
Figs. III. 4a and III. 4b respectively. Also shown by dashed lines is the 
generalization for arbitrary 6 of the "Kovner" lines given by Eq. (Ill- 14). 
These dashed lines, which correspond to the intersection points P and Q 
in Fig. Ill .2, were shown above to always lie in the region of instability. 

It is easy to determine that the minimum in the dashed line occurs at 
k C A / S2 = (4/ 3)*/^ ~ 1.15 for which 


M d 


3 (3) 1 / 2 2. 6 

2(1-5) ~ 1-5 ' 


(III- 16) 
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D 5624 


Fig. III. 3 The Dispersion Relation (III. 8) for: a) 5 = 0. 1, M d = 3; 

b) 6 = 0. 1, M d = 7; c) 6 = 0. 5, M d = 3; d) 6 = 0. 5, M d = 7. 
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M<j = U/C A 




Fig. III. 4 


Range of Unstable Values of k as a Function of Relative Drift 
Velocity for a) 6 = 0. 1 and b) 6 = 0. 5. 
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This value corresponds to the relative drift such that the beam line 
Fig. III. 2 is tangent to the W + branch, i. e. , the points P and Q coalesce. 

For larger values of drift there are two peaks in the spectrum of unstable 
k. Figure III. 4a shows that for > 3. 5, there are two discrete ranges 
of k which are unstable. For large drift (M d — *■ °° ) there are unstable modes 
for k in the range 


(1 - 6) M d 
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M d (l - 6 




kC 
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(III- 17) 


The width of this unstable spectrum is thus 


Ak 


4 6 1 / 2 
(1 - 6) U 



(III- 16) 


and the peak growth rate is given approximately by Imw ~ 6 fi These 

expressions are in agreement with the numerical results shown in Figs. III. 3 
and III. 4. 

C. RELEVANCE FOR SHOCK WAVE STRUCTURE 


The purpose of this section is to determine when the unstable whistlers 
analyzed in Section III-B can grow to large amplitude inside a shock front 
propagating along the magnetic field. Hence, in addition to being unstable 
we shall apply the requirement that the growing wave should have zero 
group velocity in the shock frame. The group velocity of the wave is not 
only a function of qj and k, but also of the position in the shock front as we 
shall see below. Waves with zero group velocity can be expected to remain 
in the unstable region for a sufficiently long time to reach large amplitude 
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and generate turbulence. We consider only parallel shock waves for which 

1 /2 

the sound speed Cg^ - (y K(T + T.)/rm) and the Alfven speed 
2 l/2 

C AO = < B o ' in the unshocked gas are equal. Thus our study 

14 

excludes any consideration of "switch-on" shocks that occur for Cg q< 
here the external magnetic field remains constant and parallel to the 
direction of plasma flow across the shock layer. 

In the rest frame of the shock front, the ions are modeled as a 
linear combination of the upstream (u) and downstream (d) velocity 
distribution functions, i. e. 


f (x, v) = n u (x) 6 (y - v ^ + m, (x) 6 (v - y H ) 


u 


(III- 18) 


njx) is taken to decrease linearly from its upstream value at the leading 
edge of the shock to zero at the trailing edge while n^(x) increases linearly 
from zero at the leading edge to its downstream value at the trailing edge as 
indicated schematically in Fig. III. 5. The electrons, on the other hand, 
are modeled as a single fluid moving with an average drift velocity, V g (x) 
determined by the zero current condition 

Ye (x) = (1 " 5(x > ) Y u + 5(x) Yd (III- 1 9 ) 

where the fraction of the total density represented by the downstream species, 
6(x), is now a function of position through the shock. 

In this section the dispersion relation (III-8) is analyzed in the shock 
frame. The primary difference from the previous section is that the real 
part of the frequency is doppler shifted by the amount k • V g (x) which is a 
function of position in the shock. The instability properties are unchanged 
and depend only on the relative drift U = (which is independent of 

position in the shock) and 6(x) (which is dependent on the position in the shock). 
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Fig. III. 5 Sketch of the Ion Velocity Distribution Function Indicating the 
Assumed Linear Variation Through the Shock. 
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The group velocity, defined as 8 Cl) / 3 k where u) is the real part of 

s s 

the frequency C0 g in the shock frame (equal to CD + k . V (x) ), also depends 
on the position within the shock. 

In the previous section, dispersion curves and stability criteria 
were given in terms of the relative drift velocity, U, of the two ion 
streams or, more exactly, the relative drift Alfven Mach number, M^ = 

i /z 

U/C^, where C A = B Q (z/ji o rm n^y' ^ is the local Alfven velocity. In this 
section we shall specify these properties and criteria in terms of the 
upstream Alfven Mach number, M^ = V /C A q, and the fraction of the 
distance, t, , through the shock. From Eq, (III- 1 8) we have 


n (x) = (1 - £ )n + £ n 
e s uo do 


(III-20) 


where n^ Q and n^ Q represent the unshocked, upstream density and 
shocked, downstream density respectively. Thus 




5 

[5 + r (1 - 5)1 


(III-21) 


where j = n ( j 0 / n u0 i s the compression ratio across the shock determined 
by the Rankin e -Hug oniot relations; 


T 


(v + 1) m a 2 

(y - l) m a ^ + 2 


4 M 2 

■= (for y = 5/3) 

m a 2 + 3 


(III -2 2) 


when it is remembered that we have assumed (3 ~ 1 so that the upstream 
acoustic speed and Alfven speed are equal. The local Alfven velocity 
decreases through the shock as the plasma density increases; 
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(III-2 3) 


C A " C AO 


M a + 3 


M a 2 + 3 + 3 UM a 2 - 1) 


]/2 


For a given upstream Mach number, M A> the local drift Mach number is 
given by 


M. 


V - 
u d 


3(M A 2 - 1) / 3^(M a 2 -1) \*/ 2 


1 + 


(III-24) 


4M, 


M a +3 


Solutions of the dispersion relation (III- 8 ) in the shock frame are 
shown in Figs. Ill .6a through III.6d for the same values of 5 and M^ 
shown in Figs. III. 3a through III. 3d. Note that the wavenumber scale 
in Figs. III. 6 are normalized using the upstream Alfven velocity C A q 
instead of the local value. The shape of the unstable spectrum is unchanged. 
The real frequencies for the case 6 = 0. 5 are no longer symmetric about 
k = 0 due to the doppler shift. Note that the group velocity does not equal 
zero (solid line with zero slope) for any of the values of k that are unstable 
in these particular cases. The case for 6 = 0. 1 and M A = 4. 1 has the 
smallest group velocity. The values of Alfven Mach number shown are 
obtained by inverting Eqs. (III-21) and (HI-24). 

Figures III. 7a and III. 7b present the range of unstable wave 
numbers as a function of the upstream Alfven Mach number, M A , at 
two different positions within the shock front. These figures are different 
from their counterparts in Figs. III. 4a and IH. 4b since the value of 5 
varies for different values of M A becuase f, the position within the shock 
is held fixed. We see that near the front of the shock there is a wide range 
of unstable wavenumbers for M A > 3. The growth rate is of the order of 
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Fig. III. 7 Range of Unstable Wavenumbers as a Function of Upstream 
Alfven Mach Number at Two Positions within the Front at 
a) £, =0. 05 and b ) =0. 2. The dashed line corresponds to 

values of k for which the group velocity of unstable waves is 
zero in the shock front. 




S2 c j. However, most of these waves have positive group velocities and 
are swept downstream with the flow before they have had a chance to grow 
significantly. Only unstable waves with nearly zero group velocity can be 
expected to grow to large amplitude. In Fig. III. 7a, the values of k for 
zero group velocity waves in the unstable region are shown by a dashed line. 
Thus, Mach numbers in the range 6. 8 < < 9. 1 should be able to produce 

turbulence at £, > 0. 05 according to the linearized analysis presented here. 
Such waves don't exist for the parameter ranges shown in Fig. III. 7b at 
£ > 0 . 2 . 

Figure III. 8 shows the range of upstream Alfven Mach numbers for 

which unstable whistlers will stand in the shock as a function of distance 

into the shock. Two immediate conclusions can be made. The first is that 

this mechanism can only be important for M^ > 5. 5. The second is that 

these waves are confined to the leading edge of the shock; unstable waves 

near the middle or trailing edge will be swept downstream. This result 

15 

is in agreement with laboratory experiments that measured magnetic 
field fluctuations in lower Mach number shock waves. 

D. SUMMARY 

We have derived the electromagnetic dispersion relation for two 
counterstreaming ion beams in a stationary electron background. The 
dispersion relation which is unaffected by electron pressure terms in 
the fluid approximation has been solved numerically to find the regions 
of instability in the k - M^ space. It is found that for beam -whistler modes 
to be unstable, it is necessary that the relative drift Mach number be 
greater than 2. 0 although this condition is not sufficient for instability. 
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Fig. III. 8 Range of Alfven Mach Numbers for Which Unstable Whistlers 
will Stand in the Shock as a Function of Position Within the 
Shock. 
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The bandwidth of k over which the unstable modes exist is shown to decrease 


inversely with the relative drift velocity for large relative drift velocities. 

By modeling the shock structure in a Mott -Smith formalism as 
consisting of streaming up-stream and down-stream ions, the linear 
analysis was applied to the study of shock wave structure propagating along 
the magnetic field. Invoking the criterion that only unstable waves with 
zero group velocity will have time to grow within the shock thickness, we 
found the range of Mach numbers for which some unstable beam-whistler 
modes could grow to large amplitudes. We suggest that this may be a 
possible mechanism to generate turbulence that many experiments have 
shown to be associated with the shock wave structures. The plasma 
turbulence in turn can provide the explanation for anomalous dissipation. 
However, in our model, the unstable beam -whistler modes can stand near 
the leading edge of the shock only for upstream Mach numbers greater 
than 5. 5. 

The most serious limitation of the present analysis is the restriction 
k| | B. For waves propagating at an angle to the magnetic field, the dispersion 
relation becomes much more complicated and electron thermal effects have 
a significant interaction with the magnetic modes treated here. The analysis 
presented here ought to be extended to include the stability properties of 
these obliquely propagating waves. The assumption of cold ions ought to be 
relaxed particularly since one of the ion species in this study has been taken 
to be the shocked ions. Ion thermal effects could have a significant effect 
on the stability properties of counter -streaming ion beams. Asa final 
comment, it should be emphasized that this instability depends on having 
included the effect of the magnetic field on the ions. If ► 0 in Eq. (III. 8), 
the instability disappears. 
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